A possible N N* (IMO) quasi- molecular state 



Lu Zhao^'^, Peng-Nian Shen^'^'"^, Ying-Jie Zhang^'^, and Bing-Song Zou^'^'^'^ 

1. Institute of High Energy Physics, CAS, Beijing 100049, China 

2. University of Chinese Academy of Sciences, Beijing 100049, China 

3. College of Physics and Technology, 

■ Guangxi Normal University, Guilin 54 1004, China 

^ ■ 

i 4- Theoretical Physics Center for Science Facilities, CAS, Beijing 100049, China 

CN 

O ' 5. Physics Department, Hebei University, Baoding 071002, China 

D : 

pLn ' 6. State Key Laboratory of Theoretical Physics, 

, Institute of Theoretical Physics, CAS, Beijing 100190, China 

' Inspired by the recent observation of a narrow resonance-like structure around 

^ • 2360 MeV in the pn — t- dir^n^ cross section, the possibihty of forming a A^A^*(1440) 
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a, p and cj exchanges in t- and u-channels are considered. By adopting the coupling 
constants extracted from the relevant NN scattering and A^*(1440) decay data, it 
is found that a deuteron-like quasi-molecular state of A^A^*(1440) with a binding 
energy in the range of 2 ~ QlMeV can be formed. Therefore, it is speculated that 
the observed structure around 2360 MeV might be or may have a large component 
of the A^Af*(1440) quasi-molecular state. 

PACS numbers: 13.75.-n, 13.75.Cs, 14.20.Gk 

I. INTRODUCTION 



The ABC effect, which shows an unexpected peak in the invariant Trvr-mass spectrum of 
the double pionic D-p fusion to ^He process, was observed by Abashian, Booth and Crowe 
long ago[l]. Later on, in the hadronic experiment, such a phenomenon was also observed in 
the other processes, for instance n + p — j- Diji-n) and D + D ^ ^HeiitTi) . The 

27r-production mechanism used to be attributed to the effect of the double A excitation. 



But it was found that the obtained enhancement in t 



TT TT channel is much smaller than the data value 



m 



le low TTTT mass region in the exclusive 



which implies that some important 
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physics is missing in such a mechanism. Thus, a sequential mechanism was proposed to 
explain the deuteron spectrum of the n + p — t- Diji-K) process, namely the nucleon can 
firstly be excited into Roper resonance iV*(1440), and then the resonance decays into iVyrTr 
and Avr subsequently [8|. 

On the other hand, the deuteron is a unique B=2 state which has been confirmed in the 
experiment at this moment, where B stands for the baryon number. Actually, after three 
decades of experimental efforts, we have never found any decisive sign for the existence of a 
B=2 state other than the deuteron. Recently, the CELSIUS-WASA Collaboration observed 
a narrow resonance^ike structure around 2360 MeV in the data set of the pn — )■ dn^ir^ 
total cross section [9]. Since such an energy is fairly close to the A^A^*(1440) threshold, 
and iV*(1440) has the same quantum numbers as the Na system and would couple to Na 
strongly, this phenomenon awakens our interest in exploring the possible existence of a 
A^A^*(1440) quasi-molecular state, which is an analogue of the deuteron. 

In order to carry out such a study, the hadron-hadron interaction should be known 
before hand. The hadron-hadron interaction used to be studied in the hadron degrees of 
freedom. Since the QCD theory emerged in the 1970s, such an interaction has further been 
investigated in the quark degrees of freedom, where the hadron-hadron system is considered 
as a multi-quark cluster system. One of the models in this category, which can basically 
explain the experiment data, is called constituent quark model (CQM). In this model, the 
quark-quark interaction can be described by the one-gluon-exchange (OGE) potential in the 
short-range, the scalar- and pseudoscalar-meson-exchange potentials in the medium-range 
interaction, and the phenomenological confinement potential in the long-range. 

However, because of the property of the non-Abelian group, the non-perturbative effect of 
QCD still cannot be taken into account accurately up to now. For simplicity and efficiency, 
the method in the framework of the hadron degrees of freedom is still widely used in the 
hadron physics investigation. In this framework, the hadron is assumed as a point-like 
particle, and the hadron-hadron interaction is supposed to be originated from exchanges 
of various mesons whose masses and coupling forms with various hadrons are different. A 
typical theoretical model is the Bonn meson-exchange model which has intensively been 
studied over years. By employing this model, a large amount of nucleon-nucleon interaction 
data, such as the deuteron property, scattering phase shifts and reaction cross sections, 
can be well described 10|]. This model has also been adopted to study the baryon-baryon 
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interaction in the coupled channel approach, although the result was not accurate enough due 



to too much approximation 



Recently, a simplified one-pion-exchange model was used 



to study systematically the deuteron-like meson-meson bound state, called deuson, some 
possible S-wave molecular-like states in the meson-meson interaction were predicted |12l . I13| . 
The interaction of A (1405) with a nucleon was studied from the viewpoint of chiral dynamics 
and a possible quasi-bound state was found jl^ . 

In this paper, the Bonn meson-exchange model is extended to the A^A^*(1440) system 
so that the binding property of the system can be studied. The paper is organized in 
the following way: The derivation of the non-relativistic A^-A^*(1440)potential is briefiy 
introduced in section II, and the model parameters are extracted from the available data 
in section III. In section IV, the numerical result and discussion are given. Finally, the 
summary is presented in section V. 

II. BRIEF FORMULISM 

The basic Feynman diagrams of the A^-A^*(1440) interaction are shown in Fig{Tl In the 
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FIG. 1: The t-channel (left) and u-channel (right) Feynman diagrams of the A^-A^*(1440) inter- 
action. 



Bonn meson-exchange model, there are six mesons ( tt, a, p, w, 5 and if) to be considered. 
Recently, by fitting the experimental data, one showed that only four of them, namely vr, a, 
p, and cj, have relatively large coupling with A^A^*(1440) [15]. Clearly, the vr- and a-meson 
exchanges would be responsible the the long-range (r^2 fm) and the intermediate-range (1 
fm^r^2 fm) interactions, respectively, and the p- and cj-meson exchanges would provide 
the core-range (r^l fm) interaction. 
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The effective Lagrangians of the meson- A^-A^* (1440) interaction for various mesons read 

jC^NN* = -ignNN*UNl5UN*r ■ TT + h.C, (l) 

^aNN* = g<TNN*UNUN*Cr + h.C, (2) 

f 

i^pNN* = -gpNN'UNl^UN'T ■ p^-- UNCr'"'UN-f ■ {df,p„ - d^p^) + h.C, (3) 

f NN* 

C-uiNN' = -gwNN*UNl^UN'UJp-- '^^—7, un(t'"' un* {d^u;^, - d^u ^) + h.C, (4) 

The forms of the effective Lagrangians of the meson- A^-A^ or meson- A^*-A^* interaction are 
the same as those of the meson-A^-A^*(1440) interaction, except the spinor. 

Now, we consider a A^A^*(1440) system whose quantum numbers are the same as those 
of the deuteron, namely I = 0, = 1+. Since A^* is also an isospin doublet, one can write 
the |A^A^*) or |A^*A^) state as 

liViv*) = ^}^\np*) + ^\pn*), (5) 

\N*N) = -^\n*p) + ^\p*n), (6) 

and its overlap as 

{NN*\NN*) = -{np*\np*) + -{pn*\pn*) — -{pn*\np*) — -{np*\pn*) , 
{NN*\N*N) = -{np*\n*p) — - {np*\p*n) — -{pn*\n*p) + - {pn*\p*n), 

{N*N\N*N) = -{n*p\n*p) + -{p*n\p*n) — -{p*n\n*p) — -{n*p\p*n) , 

1111 

(A^*A^|A^A^*) = -{n*p\np*) - -{n*p\pn*) - -{p*n\np*) + -{p*n\pn*). (7) 

Consequently, one can further write the wave function of the A^A^*(1440) system as 

1 

and the normalization as 



^) = ^( |A^A^*) + |A^*A^) ), (8) 



(^1^) = 1[ {NN*\NN*) + {NN*\N*N) + {N*N\N*N) + {N*N\NN*) ]. (9) 

Apparently, each term in the above equation contains four non-equivalent Feynman dia- 
grams. All together, Eq.(l9]) contains eight non-equivalent Feynman diagrams, as shown 
in Figj2l Among them, the first four are w-channel diagrams and the rest are t-channel 
diagrams. 
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FIG. 2: 8 non-equivalence Feynman diagrams in 

Next, we derive the N-N* potential for the above mentioned diagrams according to the 
quantum scattering theory. The relativistic S-matrix for a scattering process reads 

{f\S\t) = + i2n)'5\pf - p,)tMf,, (10) 

where Mfi denotes the scattering amplitude from the i state to the / state, or the Lorentz- 
invariant matrix element between the i state and the / state. On the other hand, by applying 
Bonn approximation on the Lippmann-Schwinger equation, the S-matrix reads 

{f\S\t) = 6f, - 27c6{Ef - E,)tVf, (11) 

with Vfi being the interaction potential. Considering the different normalization conventions 
used for the scattering amplitude Mfi and for the T-matrix Tfi, and consequently Vfi, we 
have 

Vf. = , ^ -Mf,, (12) 

f ^ i 
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where P/(i) denotes the four momentum of the final (initial) state. 

Let us use, in the center mass system, Pi{Ei,p) and P2{E2,—'p) to represent the four 
momenta of the initial particles, P^{E^,p') and P^i{E^, —p') to denote the four momenta of 
the final particles, respectively. And then, 

q^P^-P^^{Ez- El J -f) = {E2- Ei, q) (13) 

should be the transferred four momentum or the four momentum of the meson propagator. 
For convenience, we always use 

q=p'—p (14) 

and 

k^lip'+p) (15) 

instead of p' and p in the practical calculation. 
In terms of the Feynman rules, the Dirac spinor 

/ 



. , E + M 



1 

a-q/{E + M) 



6, (16) 



^ , ^+ lE + M ( \ 
^^^''^^^^UTV -^■^/{E-rM)) (17) 



and the nonrealistic reduction 



E 



the scattering amplitude caused by the 7r-meson exchange can be written as 

■M --Xn n 1 {<Ji ■ q){<J2 ■ q) 2 {mi + 1712? {(Ji ■ q){a2 ■ o) . . 

4mim2 + 8m^m2 q — 

and the scattering amplitude induced by the cr-meson exchange can be expressed by 

= -gaNN*gaN*N*^, T 1 ,22 ^ + 2 2 ^ ,22 ^'^ " X k)\ 

q^ + m% imf'm2 l&mim2 4mjm2 

- + g;;^^ - ^iS-il- k)Y (20) 

The scattering amplitudes arisen from the ui- and p-meson exchanges can be derived in the 
same way. Their forms are quite similar to those of the tt- and u-meson exchanges. To 

reduce the length of the paper, they will not be presented here. It should be mentioned 
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that because of the mass difference between and A^*, there would be a so-called " effective 
mass" Mjr,o-,p,aj in the amplitude formula 

= (m7v.(i44o) - itinY - 1^1, (21) 

^a,p,<^ = ^l,p,u^ - ("^iv*(i440) - rriNy. (22) 

Moreover, in the derivation, the term with symmetric combination (al + d^) which leads to 
a spin-orbital term S ■ {q x k) is remained, while the term with asymmetric combination 
{(fi — (J2) is dropped. Substituting M/j into Eq.(|T2l). a meson-exchange potential in the 
momentum space, V{q, k), can be obtained. 

Further making Fourier-transformation, the A^-A^*(1440) potential in the coordinate 
space, Vm(?"), can be derived. It should be pointed out that because the mass of the vr- 
meson is relatively small, the Fourier transformation of the 7r-meson-exchange potential 
would be a complex function, which implies that the off-shell effect of the exchanged pion 
would induce other reaction channels. Since what we are interested in is the elastic scatter- 
ing, as an approximation, the imaginary part of the potential can temporarily be dropped, 
which is consistent with the treatment in Ref . jl2| . For instance, 

with 

Z'(Mr) = ^^^^^^"f _ 3 sin(Mr) _ 3 cos(Mr) 

^ ^ Mr Mr Mr {Mrf Mr ' ^ ' 

Considering the finite size effect of (A^*), similar to the derivation of the Bonn potential, 
a form factor 

^, , K^ — rri^ .k^ — rri^. 

ni) - ^ - (^) (25) 

with A being the cutoff parameter and 

= A2 + m2 - (26) 

is adopted to suppress the high momentum contribution in the relatively lower energy scat- 
tering. Then the ultraviolet divergence in the Fourier transformation can be eliminated, and 
the potential in the coordinate space at small r is remained finite. Finally, the potential of 
the A^-A^*(1440) interaction can be written as 

V{r) = V^ + V, + Vp + V^ = Vc + Vls + ^r, (27) 



8 



where Vc, Vls and Vr denote the total central, spin-orbital and tensor potentials, 

. . / X gTrNNQirN^N* 3 ^ , c TT 1 

K(r)= ^^[(a,-a,)^3a + ^i2-F3.,] 

(ol ■ a.)^a.i + ^i.-Fe„. ], (28) 



9aNNgaN*N*r -rr . mj + ml ml + ml mj + m^^ - 

V'o-l.n — L ~ U "I 1 2 2~-^ 4t 77 2 2 •'St "I 2 2"'"' ' ^''5t0 

4:71 Amim2 IGmfma 4m^m2 



[ ~ J'lu + 7^ iu ~ a ^ 2« + 7^ S • LjF^uQ ], (29) 



47r 2mim2 8mim2 2mim2 



Y _ 9pNNgpN*N' r 3mf + 3ml + 12mim2 ^ ^ 3mj + 3ml j^ 

47r " im\ml Idmfml ^* 

3m^ + + 12mim2 - - 3^1-^2.^ ^ 3 



Amfml 2mim2 2mim2 

fpNN* QpNN*^ 3 6(mi+m2)=? 3(Ti-cf2,T- -rr . 35'i2_ 1 
+ 7; J'2t -J • -t^-^StO H l-^2t — -^3*1 j -^3*2 

2m\ + 2m2 47r mi m\m2 m2 m2 

, /piv*Ar* fi-piViVr 3 6(mi + m2)^ 3cfi-cf2/T- x- ^ 35'i2 ^ , 

+ -Z rr, A -^2* -J • ^-r 5to H [J- 2t — -r Ml) -''3*2 

2mi + 2777.2 47r m2 mi?772 mi mi 

9pN*Nr OX- 3"^i + 3"^2 + 12mim2 _ 3m^ + 3m2 
47r 4m(m2 l'om{m2 

3ml + 3"^i2 + 3mim2 ^ 3(mi - m2)^ ^ ^ 

5 2 ■ ^'^5u0 Z 2 2 '^l ■ 

2mfm2 Amfm2 
3{mi + m2f ^ ^, , 3(mi+m2)2 

"I T7 2 2 ■ '^2(,-^2u — •/ 3ulJ 77 2 2 '-'12-/ 3m2 

16mjm2 16mfm2 
i/piv^iv (mi - m2)2 . _3_7r 3 3 g 



An m^ 2mim2 8mim2 2mim2 



5tiO 



^ /pAfAT* fi'pAf jv* r _ 3(mi - m2f{mi + ^2) ^ _^ 3(mi + m2f ^ 
2mi + 2m2 47r mfml Amlml 

6(mi + m2)(mi + m2) ^ - 3(mi + m2)(mi - m2)^ ^ ^ 

-O • LJ- 5 o Cl ■ <^2J' Au 



mim2 mfm? 
3(mi + m2)^ ^ ^, T- \ 3(mi + m2) 



2 2 ' (^2{^2u~ ^ 3«l) ;j 2 2 *S'i2J^ 3«2 



4m^m| " " 4mf 



ms 



/piv*Ar r _ 3(mi - m2)^ _ 3(mi + m2f{mi - m2Y _ 

4w ~r ^99 *' 2u 



47r(2mi + 2m2)^ mfm^ Amlml 

6(mi + m2)^(mi - m2)^ - - 3(mi + m2)^(mi - m2)^ ^ ^ 

-O • LJ-^uQ 5 5 (Ti • (T2J-AU 



mfmg mfm2 
3(mi + m2)^ ^_ 3(mi + m2)^ ^ 

+ 4mfml '"^ ■ '"^^•^^^ " -^^"^^ " 4m?ml '^^^'^^"^ ^ 
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and 

V^(r) = —-Vp{r){ rrip — )■ m^, fpNN* — ^ JluNN*, fpN*N* — ^ fwN*N*, fpNN — ^ JluNN, 

QpNN* gujNN", gpN'N* duiN'-N", dpNN ~> QujNn) (31) 

represent the potentials induced by the vr-, a-, p- and w-meson exchanges, respectively. The 
explicit forms of J^^t^ ^au-, ^atp and J^auis are shown in the Appendix. 

As to the cutoff parameters Aq, (i.e. A in Eqs.f l25|l26|) ). it is reasonable to take the 
values used in the Bonn A^-A^ potential temporarily at the very beginning (A^r = l.^GeV , 



Ao- = l.lGeV^, Ap = l.?)GeV and A^^ = l.bGeV , which is called Set Aa) [10|, because 
A^*(1440) is an excited state of A^ and A^A^*(1440) has the same quantum numbers as the 
deuteron. However, comparing with A^, A^*(1440) has one p-wave excitation, we should 
check the rationality of this adoption, re-adjust Aq,s to some reasonable values (we will call 
it as Set B^) and use them in the later calculation. 

III. COUPLING CONSTANTS 

Before realistic calculation, 18 coupling constants should be fixed, among them gnNN, 

gaNN, gpNN, gcuNN, fpNN and f^NN for meSOU-A^-A^ couplings, fi',rAriV*(1440), fi'<7iVAf(1440), 

gpNN*(iM0), gujNN* (1440), fpNN* (1440) and /,^ArAr*(i44o) for mesou-A^-A^* (1440) couplings and 
gnN*N*, gaN*N*, gpN*N*, gujN*N*, fpN*N* and fujN*N* foT mcsou- A^* ( 1440)-A^* (1440) cou- 
plings, respectively. 

The first 6 coupling constants, g^^NN, gaNN, gpNN, g^iNN, fpNN and /ojnn, can be taken 



from the Bonn potential 10|, because the 6- and 77- meson exchanges do not contribute due 
to their strange constituents. 

The coupling constants for the meson- A^-A^* (1440) interaction can phenomenologically 
be extracted from the relevant data fitting. More specifically, (^tt at at* (1440) can be obtained 
by fitting the partial decay width of A^*(1440) to A^vr, namely the total width rAr*(i44q^ of 



200 ~ 450(^ 300)MeV times the branching ratio 5i?(A^*(1440) ^ Nn) of 0.55 ~ 0.75 ll6|. 
The corresponding Feynman diagram is shown in the left diagram in FiglSl Based on this 
diagram, the coupling constant can be derived as 



5'7r7VAr*(1440) _ 4:m^SN*^Nn /„r,\ 

An 3[(mAr* -|- rriNy — ml]^/'^[{mN' — rnNY " ml]^/'^ ' 
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FIG. 3: Feynman diagrams for the A^*(1440) A^vr (left) and iV*(1440) ^ Na ^ Ntttt (right). 



with ^Ar._^.Ar^ = rAr*(i44o) X B R{N* {IMO) — )■ Nn) being the partial decay width. Substituting 
the r7v*^Af7r data into Eg. (1521) . the value of 5'7rArAr*(i44o) can be extracted. They are tabulated 
in Table HI 



TABLE I: gaNN*{i440)^ faNN*(i440) ^^'^ with a = iT{138MeV), a{550MeV), p{775MeV) and 
u;{782MeV) in Set ^a- 



Coupling mode of p 


Coupling constant 






a 




vr 


a 


P 




V+T 


daNN* (1440) /^^ 


2.9087 


5.1200 


0.3685 


1.3890 


faNN* (1440)/^'^ 






2.8288 




V 


9 aNN* (1440)/'^'^ 


2.9087 


5.1200 


1.1374 


1.3890 


faNN* (1440)/'^'^ 










T 


9aNN* (1440)/^'^ 


2.9087 


5.1200 




1.3890 


faNN*{1440)/^'^ 






1.1374 




V+T 


alNNl^-^ 


14.900 


7.7823 


0.9500 


20.000 








35.350 






K (GeV) 


1.3000 


1.1000 


1.3000 


1.5000 



In the same way, 5'o-ArAr*(i44o) can be extracted by fitting the branching ratio 
i?i?(iV* (1440) — )■ Na — )■ N-mx) of 0.11 ~ 0.13 [16]. The corresponding Feynman diagram is 
shown in the right diagram in FigJHJ Based on this diagram, the coupling constant can be 
written as 

9aNN*{l440) _ ^ N'^Ntttt Sn^misf* , . 
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with factor S being 

^= / 72 1 ^ Y2T( ^ ^2 2\2^ 2^ g Q ^ [P ^ Pi ^ PZ - Pa) ■ (34) 

J - [Ps + PiY {{P3 + PiY - mlY + mlVl ^p\plpl 
Extracting the value of the couphng constant for the p-meson is somewhat comphcated. 
There are two types of couphng modes in the p-N-N* interaction, thus two couphng con- 
stants, the vector type couphng constant gpNN*{i4:40) and the tensor type couphng constant 
/piv Af* (1440)) should be fixed. To do so, we need two experimental values. Here, we choose 
the partial widths of the A^*(1440) Np Ntitt and A^*(1440) A^7 (/^ = 1) decays. 
The former width can be obtained by multiplying the total width Tm^i^^q-j of about 0.3 



GeV and the branching ratio BR{N*{U40) Np ^ Ntttt) of < 8% [16|. In the practical 
calculation, taking BR{N* {IMO) Np Nnn) = 2% would be reasonable. The later 
width can be calculated by the following equation 



k'^ rri]^ 



Tn-^n^ {i,=i) = ■^—r(^i]2 - KI2) (35) 

with A^j^ = -0.065 GeV-^/^ and A'lJ^ = 0.04 GeV~^/^ from the PDG Q, respectively. 
The corresponding Feynman diagrams are shown in FigHl where the vector-meson-dominant 
mechanism is employed in the right diagram. Based on this figure and FejTiman rules, we 



N 

N*{UAO) 




N 

N*{U40) 




FIG. 4: Feynman diagrams for the iV*(1440) ^ Np^ Ntttt (left) and iV*(1440) p-y {I^ = 1) 
(right). 



can write the partial decay widths for these processes as 
Ftv-j-TYtt+tt- = / 5^{P — Pi — P3 — Pi 



4.0 o .d^Pi d^p3 d^pi TUN A -ml 



2P0 2P3O 2P0 (27r)5 ' A2 - p2 / (p2 _ ^2)2 

fpNN 

rriN* +m]\r 

and 



glAiOpNN' + fpNN^fc, + (2-^^^^P ■ (P4 - P3)fc2] (36) 



i Af.^TV^ (7^=1) - — ^^^^3- { ) ^{gpNN- + JpNN*) {A ( ) 
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where the factors c,- can be written as 



\2 



and 



c, = ^^^^ + 1, (39) 

■m]\i*m]\i 



respectively, and the p7 couphng constants gp^ is given by [17 1 



with fp/4:7i = 2.7. Solving coupled Eqs. (l36|l and (1371) . we can obtain the values of gpNN*{iMo) 
and /pArAr.(i44o). We tabulate them in Table [T] as well. 

To get more information about the p-N-N* coupling, we also arrange for the cases in 
which only one p-N-N* coupling mode, either the vector or the tensor, exists. In these 
cases, one value of the partial decay width would be enough to fix either 5'pArAr*(i44o) or 
fpNN*{u40)- Then, the corresponding partial width formulas become 



N-'-^N-^l (1^ = 1) - QpNN' 



47r 2m^. ' mf 



(41) 



and 



P _ ( fpNN* ^23% {mlf, - m%) {niN* - mjy)^ ^ A^ - 1 

i N'^N, (/,=i) - l2m^.+2m^^ 4^ ^^V^^ 



(42) 



Substituting the measured partial decay width into the above equations, we obtain the value 
of fl'pAfAr*{i440) oi fpNN*(iAm)- We also list them in Table [H 

Finally, we extract the coupling constant of the (X'-A^-A^*(1440) interaction. As well 
known, the u-N-N interaction has a vector type coupling only. Due to lack of Ftv-s-atcj 
data at the present moment, roughly taking a single vector coupling mode for the u-N- 
A^*(1440) interaction would make sense. Thus, we can use the partial decay width of the 
AT* (1440) ^ N-f (/^ = 0) process with A^]^ = -0.065 GeV^^^^ and A^J^ = 0.04 GeV-^/^ Q 
to determine gujNN*{i4m)- The corresponding Feynman diagram is shown in Figj5l According 
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A^*(1440) 




FIG. 5: Feynman diagrams for the Af*(1440) — > {I^ = 0). 



to this diagram, we can calculate the partial decay width by using 

r^..^.a.^o)= ^^iAl], + A-],r. (43) 



Meanwhile, we can write the partial decay width formula as 

r^.^^, (,^=0) - (^^) — (44) 



with 



2 eV47r • 

9.. = ^^. (45) 



In the last equation, we roughly take f'^j^T^ = df^/An = 24.3 via SU(3) relation. By fitting 
measured partial decay width, we can extract the value of 5'aj7VAf*(i44o) approximately, and 
also tabulate it in Table HI 

However, we do not have any experimental data at all to determine the coupling constants 
of the meson- A^*-A^* interaction at the present stage. Therefore, we have to surmise a way 
to set the values of gaN*N*, (a = vr, cr, p, u) and fi3N*N*, (/3 = P, w). Here, we consider 
two cases: Case 1, gaN*N* = 9aNN, (a = vr, a, p, u) and ff}N*N* = ffiNN, (/? = P, uj); and 
Case 2, gaN*N' = 9aNN', (tt = TT, a, p, u) and fpN'N' = fpNN*, {P = P, w)- In each case, 
we have three coupling modes for the p — N — N* interaction, namely the vector+tensor 
coupling mode(VT), the vector coupling mode only (V) and the tensor coupling mode only 
(T). 

Based on the uncertainties of the data, both the total width of A^*(1440) and the mea- 
sured branching ratios for the above mentioned processes, we further tabulate corresponding 
partial widths and, consequently, the extended ranges of the extracted coupling constants 
in Table M 
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TABLE II: Uncertainties of partial decay widths in concerned decay processes and corresponding 
variational ranges of gaNN*{i4:40) faNN* (1440) ^or vr, a, p and uj with the values of in Set 



a 


TT 


(J 


P 




branching ratio 


r./r 




^N'^N-y {Ij=l) 




0.55 ~ 0.75 


0.11 ~ 0.13 


0.074 ~ 0.13 


0.0011 ~ 0.014 


dltMN* {1440)1^'^ 


2.461 ~ 3.356 


4.693 ~ 5.547 


0.8543 ~ 1.461 


0.2689 ~ 3.380 


faNN'-{144{))l^'^ 






0.8543 ~ 1.461 





To get a comprehensive impression, we plot the upper and lower bounds of the potentials, 
caused by the uncertainty of the partial decay width data, in the cases IV and 2V in Figj6l 
in the cases IT and 2T in FigJTJ and in the cases IVT and 2VT in FigJHl respectively. 
The solid and dashed curves denote the upper and lower bounds of the central potential, 
while the dotted and dash-dotted curves represent the upper and lower bounds of the tensor 
potential. From these curves, one sees that for the central potential, although a strong 
repulsion appears in the small r region, an even stronger attractive effect in the medium 
and large r regions might bring and A^*(1440) together to form a bound state. The 
tensor potentials in the case 2, no matter with which p coupling mode, are stronger than 
those in the case 1, thus the D-state might couple to the S'-state strongly, which is also 
beneficial for forming a bound state like in the deuteron case. Moreover, the strengths of 
the potentials in the cases with the T and VT coupling modes are closer but differ with that 
in the case with the V coupling mode, which implies that the tensor coupling mode in the 
p-A^-A^*(1440) interaction plays a more important role than the vector mode does. Thus, 
the binding properties of the A^-A^*(1440) system in the T and VT cases would be closer 
and would deviate from those in the V case. 

We also plot the upper and lower bounds of the central potentials contributed by various 
mesons in the case IV in Fig. O The solid, dashed, dash-dotted and dotted curves denote the 
contributions from the tt-, a-, p- and cu-mesons, respectively. Clearly, the a-meson provides a 
major attractive force which plays a key role in binding, the p-meson gives a weak attraction 
in the medium r region and a strong repulsion in the small r region, the contribution from 
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FIG. 6: iV-A^*(1440) potentials in the cases IV (left) and 2V (right). The solid and dashed curves 
denote the upper and lower bounds of the central potential, while the dotted and dash-dotted 
curves represent the upper and lower bounds of the tensor potential. 




FIG. 7: A^-A^*(1440) potentials in the cases IT (left) and 2T (right). The solid and dashed curves 
denote the upper and lower bounds of the central potential, while the dotted and dash-dotted 
curves represent the upper and lower bounds of the tensor potential. 



TT-meson shows a very weak attractive feature in the long distance and a strong repulsive 
character in the short and medium distance, and the cj-meson mainly provides a repulsion. 
The corresponding characters in the T and VT cases are similar to these, we would not show 
them again. 
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FIG. 8: iV-iV*(1440) potentials in the cases IVT (left) and 2VT (right). The solid and dashed 
curves denote the upper and lower bounds of the central potential, while the dotted and dash-dotted 
curves represent the upper and lower bounds of the tensor potential. 




FIG. 9: The upper and lower bounds of the central potentials contributed by various mesons for 
the case IV. The solid, dashed, dash-dotted and dotted curves denote the contributions from the 
IT-, a-, p- and cj-mesons, respectively. 



IV. RESULT AND DISCUSSION 

Having the A^-A^*(1440) potential V{r) in Eq.( l27|) . we are able to study the binding 
property of the A^-A^*(1440) system by solving the Schrodinger Equation 

(-—V' + Vir) - E)^<{r) = 0, (46) 
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where \I/(r) is the total wave function of the system. The quantum numbers of the system 
we are going to study are the total spin S = 1 and the orbital angular momenta L = and 
L = 2. Thus the wave function \I^(r) should have a form of 

"ifif) =i^s{r)+i^D{r), (47) 

where V's'(^ ^ind ipni'^ are the S'-wave and D-wave functions, respectively. In the matrix 
method, we use Laguerre polynomials as a set of orthogonal basis 



Xni[r) 



r(2/ + 3 + n 
with a normalization condition of 



poo 

/ Xim{r)Xin{r)r'^dr = SijSmn, (49) 
and expand the total wave function as 

n— 1 n— 1 

^{r) = J2 aiX^o{r)<|)s + bpXp2{r)(pD, (50) 

1=0 p=0 

where (ps and (pD are the angular part of the orbital wave function and spin wave function 
for the S- and D-states, respectively, and and 6j are corresponding expansion coefficients. 

In the practical calculation, Eq. fH6|) is re- written, by detaching the terms related to the 
kinetic-energy-operator (V^) from V{f), as 

(-— - —[V^a{r) + a(r)V^] + V{f) - E)^(f) = (51) 

^ Li ZiLL 

with _^ 
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a\r) — \—Z^)\\ 2 2~*^4i2 H -r ^u2 

An Amfm2 An 2mim2 

9pNN9pN'N* 3m^ + 3ml + 12mim2 



47r 4mf 



4t2 



91n'N , Smf + 3ml + 12mim2 3(mi - m2)^ ^ ^, 
47r Am\ml Amfml ^ ^ "^"^ 

5'pAr*jv (mi - ma)^ 3 ^ 
An 2mim2 

/pjVAf* 3(mi - m2Y{mi + m2) ,^ ^ ^ 

2mi + 2m2 An m\ml ^ ^ "^"^ 

/pAf*Ar ^3(mi-m2)'' , 3(mi + m2)^(mi - m2)^ 



47r(2mi + 2m2)^ ^ mfml ~^ mfml ^)-^4m2 

_lr 9uNN9icN*N* 3m^ + 3m| + 12mim2 ^ 

3 An Amfml ^^"^ 

I 91n*n ^ 3m^ + 3ml + 12mim2 _ 3(mi -7722)^ ^ _ 

47r Am\ml Am\ml ^ ^ "^"^ 

4n2 



47r 2mim2 
^iVTV* 5'a;iV7V 3(mi - m2)^(mi + 7712) 



2mi + 2m2 47r mfm| ( + 1 ^) 4?i2 

/(^AT'TV ,3(mi-m2)'^ , 3{mi + m2f{mi- m2Y 



H ■ 0^2)-^4«2 ]},(53) 



An{2mi + 2m2)^ m\ml m\ml 
and the effective potential 

\/(r) = ^(r) - {terms related to V^). (54) 
Then, with the wave function in Eq. (l50|) . the Hamiltonian matrix can be expressed as 



jjSS jjSD 
jjDS jjDD 



(55) 



with 



roo ^2 ^2 

= (05| / y2aiXio{r){- — [l + a{r)]V\xjo{r) - —V^[a{r)ajXjo{r)] 
Jo i , ^/^ 

+V^55(r)a,X,o(r)}r2dr|05), (56) 

poo 
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H'''={<Pd\ y2^Xv2{r)VDs{r)aa^o{rydr\ct)s). (58) 

and 

/oo »-2 fc2 

E + «(r)]V^6,x.2(r) - ^V^[a(r)6,x.2(r)] 

+VDD{r)\x,2{r)}r''dr\(PD). (59) 

In the above equations, Vssir), Vsnir), VDs(r) and VoDir) are taken to be V^(r), VT(r), 
Vt(?") and Vc(r) + Vt(?"), respectively. Diagonalizing this matrix, the binding energy and 
corresponding wave function of the A^A^*(1440) state can be obtained. 

The resultant binding energies for the S'-states with the coupling constants shown in 
Table |T] (with Set Ajs) are tabulated in Table lllli From this table, it is seen that with such 



TABLE III: Binding energies of the bound S'-states of the 7VA^*(1440) system (in MeV) with Set 



CoupUng mode of p 




Binding 


; energy 




case 1 


case 2 


V 


El 


-373.7 


-583.3 


E2 


-23.2 


-68.0 


T 


El 


-316.5 


-422.5 


E2 


-11.0 


-29.0 


V+T 


El 


-324.6 


-419.2 


E2 


-13.0 


-29.5 



a potential, no matter in which case and with which coupling mode, there might exist two 
bound S'-states. One of them is a deeply bound state whose energy is close to the energy of 
the deuteron, and the other is a weakly bound state. 

Contributions from various mesons to the binding energies of the obtained bound states 
in various cases with Set are shown in Table IIVI For simplicity, relatively smaller con- 
tributions from the D-state are disregarded. From this table, it is shown that the major 
attraction for binding comes from the cr-meson exchange, and the main repulsion for unbind- 
ing comes from the cj-meson exchange and the kinetic energy. The averaged contributions 
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TABLE IV: Contributions from various mesons and kinetic energy to the binding energies of the 
obtained bound states in various cases (in MeV) with Set A/<^. 









Binding 


; energy 








CclSG 




total 


kinetic energy 


vr 


(T 


P 




o c o I \/ — 


El 


-339.8 


146.7 


53.53 


-821.4 


-31.17 


312.5 


Odot: 1 V 


E2 


-18.74 


101.4 


13.00 


-198.9 


4.975 


60.83 


r* Q G y \/ ~ 


El 


-477.3 


239.6 


83.75 


-892.2 


-55.46 


146.9 


Cdot: Z( V 


E2 


-48.50 


186.1 


14.34 


-273.1 


-7.307 


31.47 


CO n 1 ~ 


El 


-284.5 


132.1 


36.38 


-707.3 


-6.001 


260.3 




E2 


-9.076 


56.46 


6.346 


-104.8 


3.551 


29.39 




El 


-363.2 


210.3 


56.74 


-771.2 


-19.06 


121.8 


case 2T 


E2 


-22.74 


112.4 


8.051 


-165.5 


4.256 


18.04 




El 


-297.8 


131.9 


36.44 


-715.7 


-15.71 


265.3 


case IVT - 


E2 


-11.11 


65.35 


7.407 


-123.6 


3.877 


35.89 




El 


-379.9 


182.7 


45.34 


-731.2 


7.715 


115.5 


case 2VT - 


E2 


-25.44 


108.5 


7.014 


-165.0 


5.619 


18.40 



from the vr-meson exchange, which also shows the repulsive character, and from the p-meson 
exchanges, whose character depends on the case, are much smaller than those from the a- 
and cj-exchanges. 

As mentioned previously, before realistically studying the binding property of the 
iVA^*(1440) system, the cutoff values should be carefully re-adjusted due to the larger size 
of A^*. In other word, the obtained result should be stable with respect to the model pa- 
rameters, namely the deviations of the cutoff parameter values, thus the potential model 
would have prediction power. Here, we plot the binding energy E2 as the function of Aq, 
{a = cr(upper left graph), (upper right graph), vr (lower left graph) and p(lower left graph)) 
in Fig. [TOl The curves for Ei are very similar, so we will not present them here. From this 
figure, one sees again the above mentioned features: (1). The binding energy in the single 
tensor couphng case is close to that in the vector -l-tensor coupling case, which implies that 
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cutoff of o(MeV) 

1400 1600 




cutoff of (o(MeV) 

1400 1600 





FIG. 10: Binding energy E2 as a function of Aq, (a = o"(upper left graph), a;(upper right graph), 
7r(lower left graph) and /3(lower left graph)). 

we can use a single tensor coupling mode instead of the vector+tensor coupling mode for 
simplicity. (2). The binding energy in the single vector coupling case is larger than that 
in the single T or V+T cases in the Case 2. (3). The binding energy in the Case 2 is 
larger than that in the Case 1 due to more attractive feature of the potential in the Case 
2 than in the Case 1. Besides, there exist different features for various meson exchanges. 
For the a-exchange (upper left diagram in Fig. ITOl) . we have following observations: (1). 
The binding energy increases with increasing because of an attractive feature of the a- 
exchange. (2). The result in the Case 1 is stable along with the variation of A^.. But, in the 
Case 2, the result increases rapidly with increasing A^.. This is because that the averaged 
contribution from the a-exchange is much larger than those from the other exchanges, thus 
the inconsistency between the coupling constants of a symmetric N*-a-N* interaction and 
those of an asymmetric N-a-N* interaction would be enlarged. It also implies that the 
results in the Case 1 might be more meaningful. For the cj-exchange (upper right diagram 
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in Fig. [10]), it shows: (1). The binding energy decreases with increasing A^^ because of a 
repulsive feature of the cu-exchange. (2). The result is stable as long as > l.SGeV. (3). 
When Atj < 1.2GeV, the result would be unstable, because too much repulsion is cut due 
to a smaller A^ value. For the vr-exchange (lower left diagram in Fig. [TOj) . we see: (1). The 
binding energy increases with increasing A^r in the smaller A^r region. (2). The result is 
relative stable as long as A^ ^ 1.1 ~ l.SGeV. For the p-exchange (lower right graph of 
Fig. [T0|) . we get: (1). The binding energy firstly increases then decreases with increasing Ap, 
because the character of the average contribution from p-exchange depends on how much 
repulsion being cut. (2). The result is stable when Ap ~ 1.2 ~ 1.5GeV only. Form above 
analysis, we find that in order to have a meaningful result, we should take the values of A^ 
within the following ranges: A^ ~ 1.1 ~ l.SGeV, A^- > l.lGeV, Ap ~ 1.2 ~ 1.5GeV and 
A^ > l.SGeV. Therefore, in the rest of the calculation, we take A^^ = l.lGeV, A„ = 1.2GeV, 
Ap = l.SGeV and A^ = l.SGeV, respectively. With this set of cutoff parameters, called 
Set -Ba, we re-calculate all the coupling constants and tabulate them in Table IVl The corre- 
sponding binding energies are tabulated in Table IVD From this table, we again see that the 



TABLE V: gaNN*{U40), faNN*{U40) and Aa in Set B\. 



Coupling mode of p 


Coupling constant 






a 




vr 


a 


P 




V+T 


fl'aAriV*(1440)/^^ 


2.9087 


4.900 


0.3685 


1.1187 








2.8288 




V 




2.9087 


4.900 


1.1374 


1.1187 


/aAfAf(1440)/^'^ 










T 


9aNN*{lU0)/'^'^ 


2.9087 


4.900 




1.11870 


faNN'-{lUO)/'^'^ 






1.1374 




V+T 




14.900 


7.7823 


0.9500 


20.000 








35.350 






Aa (GeV) 


1.1000 


1.2000 


1.3000 


1.8000 



binding energies in the case with a single vector coupling mode are usually larger than those 
in the case with a single tensor coupling mode or a vector+tensor coupling mode, and the 
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TABLE VI: Binding energies of the A^iV*(1440) system (in MeV) in Set B^- 



Coupling mode of p 




Binding 


; energy- 




case 1 


case 2 


V 


El 


-296.5 


-596.6 


E2 


-10.2 


-68.1 


T 


El 


-253.0 


-440.1 


E2 


-4.7 


-30.0 


V+T 


El 


-261.2 


-436.4 


E2 


-5.7 


-31.2 



TABLE VII: Variational ranges of coupling constants in Set B/^. 





TT 


(7 


P 






r./r 


r./r 






Branching ratio 


0.55 - 0.75 


0.11 - 0.13 


0.074- 0.13 


0.0011 - 0.014 


5a AT AT •(1440)/^"' 


2.461 - 3.356 


4.492 - 5.308 




0.2166 - 2.722 








0.8543- 1.461 (T case) 






2.461 - 3.356 


4.492 - 5.308 


0.02282 - 0.9717(VT case) 


0.2166 - 2.722 


faNN' (1440)/'^''^ 






3.779 - 2.084(VT case) 




A„(GcV) 


1.1 


1.2 


1.3 


1.8 



binding energies in the later two cases are comparable. For self-consistency, we would ignore 
the single V coupling mode in the later calculation. Moreover, the results in the Case 2 are 
generally larger than those in the Case 1. All these observations agree with the potentials 
shown in the previous section. 

We also present the variational ranges of coupling constants with Set B\ due to the 
uncertainties of the decay data for A^*(1440) in Table IVIII and corresponding ranges of 
binding energies in Table IVIIIi From this table, one sees that because of many uncertain 
factors, for instance the large uncertainties of the A^*(1440) decay data, the indeterminate 
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TABLE VIII: Binding energies of the iViV*(1440) system (in MeV) in Set B{^. 



Coupling mode of p 






Binding 


; energy 




case 


1 


case 2 


T 


^1 


-222.3^- 


-285.0 


-32-1.2- -581.5 


E2 


-2.34~- 


-8.330 


-12.16 62.61 


V+T 


El 


-230.2~- 


-293.7 


-330.9~ -562.6 


E2 


-2.897-- 


-9.987 


-13.38 67.43 



coupling mode of the p-meson, the coupling constants of gN*N*M and fN*N*M and etc., the 
resultant binding energies spread in the large regions, say about 222 ~ 584MeV for the 
ground state and about 2 ~ 67MeV for the first excited state. 

Before closing this section, we would mention that in the N-N scattering, nucleon might 
be excited to A^*. So, the deuteron could be a mixture of the 5"— wave ground state with 
other high partial wave states, say the D-wavc state, and even A^A^* or N*N* states, as long 
as the quantum numbers are allowed. Because the deeply bound state of iViV* (1440) with a 
binding energy of about a few hundred MeV is close to the energy level of the deuteron, it 
might couple to the conventional deuteron to form a realistic deuteron. Namely, this state 
would be a component of the real deuteron. Therefore, one may not be easy to observe 
it in experiment. However, the second bound state just has a binding energy of about a 
few tens of MeV, say 2MeV ~ 67MeV in this calculation, it is possible to be observed in 
experiment. In fact, this state is very close to the newly observed resonance-like structure 
around 2360 MeV in the pn — >■ ci7r°7r° reaction cross section data. Therefore, we speculate 
that the observed state may have a A'"A'"*(1440) structure. 

V. SUMMARY 

The meson-exchange model is extended to study the A^A^*(1440) interaction. In this 
model, the tt-, a-, p- and cj-meson exchanges between N and A^*(1440) are assumed. By 
calculating t- and w-channel Feynman diagrams, the A^-A^*(1440) potential is derived. The 
coupling constants QaNN/^'^ faNN/^'^ taken from the Bonn potential model. The 
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coupling constants 9 ^nn* {1440)/^'^ faNN*{i44Q)l^'^ phenomenologically determined 
by fitting the partial decay width data of the N* Nn, N* Nnn, N* ^ Np ^ 
Ntttt and N* — >■ A^7 processes. As to the coupling constants gaN''{iM0)N*(i440)/^'^ ^^'^ 
/«Af(i440)7V*(i440)/^'^' ^^'^ cases are assumed. It is found that similar to the N-N poten- 
tial, the (j-meson exchange provides a major attraction, the cj-meson exchanges mainly sup- 
plies a repulsive force, tt- and p-meson exchanges give an effective medium- and long-range 
attractions. Solving the Schrodinger equation by a matrix method, the binding energies 
and corresponding wave functions are obtained. It is shown that there might exist two N- 
A'"*(1440) S-wave bound states. One of them is deeply bound with a binding energy ranged 
from 222MeV to 584 McV. This state could be a small component of the real dcutcron. 
Another state is weakly bound with a binding energy ranged from 2MeV to 67 McV. This 
state is very close to the newly observed resonance-like structure around 2360 MeV in the 
pn — > dn'^n^ reaction cross section. Therefore, we speculate that the observed state may 
have a large component of the A'"A^*(1440) quasi- molecular structure. 



VI. APPENDIX 



A^-m^ 1 



= mY{mr) - KY{Kr) - (A^ - w?) ■ 



-Ar 



2A 



(60) 



^lu — 



A2 

Hi-: 



m 



)- 



MY (Mr) - Ar(Ar) - (A^ 



— m^)- 



Ar 



2A 



(61) 



2t 



Hi 



Is? - w? 



)- 



^2 _|_ g2 g2 _|_ j^2 



} 



= m^[Ay(Ar)-my(mr) ] + (A^-m^)A- 



-Ar 



(62) 



= M^[AF(Ar)-MF(Mr) ] + (A^-m^)A 



1 e"^'' 
-cfi ■ (J2[ m^AY{Ar) - m^Y{mr) + (A^ - m^)A— — ] 

+-5i2[-m=^Z(mr) + A3Z(Ar) + (A' - m')(l + Ar)^y(Ar) ] 
3 2 



A^-m^ (ai-g)(a2-g) 

-ol • M2Ay(Ar) - M^Y{Mr) + (A^ - m'^)A ] 

3 2 

+-Si2[-M^Z{Mr) + A^Z(Ar) + (A^ - m''){l + Ar)^y(Ar) ] 
3 2 



= -^{l A2 I -2 2} 

A^ + g^ + 

m^,^, , A^,^,. , A^-m^,Ar .e"^*" 
.Y{mr) - —Y{Ar) (— - 1) 



4 ' ' 4 ' ' 4 ' 2 ' r 
--{V\mY{mr)-AY{Ar) —} 

-^4*1 + {-^V2, 7-4*2} 



2 _ ^2 ^-Ar 



e 



— F(Mr) - -F(Ar) - 1)— - 

1 . . A^ - e^ '^'' 
--{V^MF(Mr)-Ar(Ar) ^ ^} 
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5t 



S-L[ - m^Z^imr) + A^Z^{At) + (A^ - m^)- 



-Ar 



2r 



S-LT, 



5tO 



(68) 



^f.i ^^-m\S- jqx k) ^ 



-Ar 



= S-L[ -M^Zi{Mr)+A^Zi{Ar) + {A^ -m'')- 



2r 



5uO 



(69) 



6u 



p2 -M2 

(tI • - M^Ay (Ar) - M^^^^^ ' ^ ^ 2 -2 



A2 + g2 

I- - 
3 

1 



+ (A^ - m')A- 

A, 



-Ar 



^12 [M3^'(Mr) + A=^Z(Ar) + (A=^ - m^){l + Ar)^F(Ar) 



{(7l ■ (J'2)Tqu\ + Svi^^t 



6u2 



(70) 
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